Preprint typeset using LAT^X style emulatcapj 



MAGNETIC PROPELLER OUTFLOWS 
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ABSTRACT 

A model is developed for magnetic 'propeller'-driven outflows which cause a rapidly rotating mag- 
netized star accreting from a disk to spin-down. Energy and angular momentum lost by the star goes 
into expelling most of the accreting disk matter. The theory gives an expression for the effective Alfven 
radius Ra (where the inflowing matter is effectively stopped) which depends on the mass accretion rate, 
the star's mass and magnetic moment, and the star's rotation rate. The model points to a mechanism 
for 'jumps' between spin-down and spin-up evolution and for the reverse transition, which are changes 
between two possible equilibrium configurations of the system. In for example the transistion from spin- 
down to spin-up states the Alfven radius Ra decreases from a value larger than the corotation radius to 
one which is smaller. In this transistion the 'propeller' goes from being "on" to being "off." The ratio of 
the spin-down to spin-up torque (or the ratio for the reverse change) in a jump is shown to be of order 
unity. 

Subject headings: accretion, accretion disks — plasmas — magnetic fields — stars: magnetic fields — X-rays: 
stars 



1. INTRODUCTION 

Observations of some X-ray pulsars show remarkable 
'jumps' between states where the pulsar is spin-ning-down 
to one where it is spinning-up. Examples include the ob- 
jects Cen X-3 (Chakrabarty et al. 1993) and GX 1+4 
(Chakrabarty et al. 1997; Cui 1997). The theoretical 
problem of disk accretion to a rotating magnetized star 
has been discussed in many works over a long period 
(Pringle & Rees 1972; Lynden-Bell & Pringle 1974; Ghosh 
& Lamb 1979; Wang 1979; Lipunov 1993; Shu et al. 1994; 
Lovelace, Romanova, & Bisnovatyi-Kogan 1995 (hereafter 
LRBK); Li & Wickramasinghe 1997). However, except for 
the work by Li & Wickramasinghe (1997), the studies do 
not specifically address the 'propeller' regime (Illarionov 
& Sunyaev 1975) where the rapid rotation of the star's 
magnctosphere acts to expell most of the accreting matter 
and where the star spins-down. Recent computer simula- 
tion studies of disk accretion to a rotating star with an 
aligned dipole magnetic field (Hayashi, Shibata, & Mat- 
sumoto 1996; Goodson, Winglee, & Bohm 1997; Miller 
& Stone 1997) provide evidence of time-dependent out- 
flows but do not give definite evidence for a 'propeller' 
regime with spin-down of the star. The present work con- 
siders the 'propeller' regime and develops a simple physical 
model where the energy and angular momentum lost by 
the rotating star goes into a magnetically driven outflow. 



2. THEORY 

We consider the problem of disk accretion onto a rotat- 
ing magnetized star which has an aligned dipole magnetic 



field. We focus on the limit where the star is rotating 
rapidly and the disk-star configuration is as sketched in 
Figure 1 (see Figure 3 of LRBK) . We consider the flow of 
mass, angular momentum, and energy into and out of the 
annular region indicated by the box A'ACC in this figure, 
where A' A is at radius r\ and CC is at r 2 . Notice that at 
this point the values of r\ and r 2 are unknown. They are 
determined by the physical considerations discussed here. 

Consider first the outer surface CC through the disk. 
The influx of mass into the considered region is 



M„ 



27rr 2 / dz [p (-v r )} 2 > 



(1) 



where h is the half-thickness of the disk, and the 2 sub- 
script indicates evaluation at r = ri. We assume that mass 
accretion rate for r > r 2 is approximately constant equal 
to M accr . That is, we consider that outflow from the disk 
is negligible for r > r<i. The influx of angular momentum 
into the considered region is 
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where T£™ is the viscous contribution to the stress ten- 
sor which includes both the turbulent hydrodynamic and 
turbulent magnetic stresses. The influx of energy into the 
considered region is 



CM 
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(3) 
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Fig. 1. — Geometry of disk accretion to a rapidly rotating star with an aligned dipole magnetic field. Here, r\ and r% are 
the boundaries of the region considered in the text; Ra is the effective radius where the outflow begins; w„ is the star's 
rotation rate; ujk is the Keplerian rotation rate of the accretion disk; and r cr = (GM /w^) 1 / 3 is the corotation radius. The 
magnetic field in the vicinity of r\ has an essential time-dependence owing to the continual processes of stellar flux leaking 
outward into the disk, the resulting field loops being inflated by the differential rotation (LRBK), and the reconnection 
between the open disk field and the closed stellar field loops. The dashed lines marked by the letters 'n' indicate neutral 
surfaces along which reconnection occurs. 



c s <^v K , where c s is the 



where u> = v^/r w vk/t with vk = \J GM/ r the Keple- 
rian speed and M the mass of the star, and where wi is 
the enthalpy. For conditions of interest here the disk at r 2 
is geometrically thin so that w 
sound speed. 

Consider next the fluxes of mass, angular momentum, 
and energy across the surface AA' in Figure f. For the 
physical regime considered, where VK(ri)/ri <C (= 
the angular rotation rate of the star), the mass accretion 
across the AA' surface is assumed to be small compared 
with M accr . The reason for this is that any plasma which 
crosses the AA' surface will be 'spun-up' to an angular 
velocity u>* (by the magnetic force) which is substantially 
larger than the Keplerian value, and thus it will be thrown 
outwards. Thus the efflux of angular momentum across 
this surface from the considered region is L± = —T% . The 
efflux of energy across this surface is E\ = —lu*Ti, where 
tt>* is the angular rotation rate of the star and the inner 
magnetosphere as shown in Figure f . For the conditions 
considered, the star slows down and loses rotational energy 
so that Ti > 0. 

We have E\/L\ = cj*. Because the interaction of the 
star with the accretion flow is by assumption entirely 
across the surface A A' ', this is consistent with the spin- 
down of a star with constant moment of inertia /; that is, 
AE^/AL* = J,w»Aw,/(/»Aw f ) = cj*. 

Next we consider the mass, angular momentum, and en- 
ergy fluxes across the surfaces A'C and AC in Figure f . 



As mentioned, accretion to the star is small for r cr <C r± 
where r cr is the corotation radius as indicated in Figure 
f. Thus, the mass accretion goes mainly into outflows, 



Mr, 



M n 



where "out" stands for outflows. The an- 



gular momentum outflow across the surfaces A'C and AC, 
L out , must be the difference between the angular momen- 
tum lost by the star and the incoming angular momentum 
of the accretion flow. The angular momentum carried by 
radiation from the disk is negligible because (vk/c) 2 <C 1, 
where c is the speed of light. That is, L out = L accr — L\. 
The energy outflow across the A'C and AC surfaces is 
E ut + E rad = E accr - Ei, where E rad is the radiation 
energy loss rate from the disk surfaces between r = r\ 
and 7"2, and E out is the rate of energy loss carried by the 
outflows. 

Angular momentum conservation gives 



Lout = M accr (rv rj> )2 -T 2 +Ti 



(4) 



We have 



• • GM Mgccr „ , „ 

J^out + J^rad = ~ ^2 -t 2 + 1 1 . (b ) 

We can solve equation (4) for T2 and thereby eliminate 
this quantity from the energy equation (5). This gives 

• • 3GM Mgccr , f , ( . T , fi , 

2r 2 
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The preceeding equations are independent of the na- 
ture of the outflows from the disk. At this point we con- 
sider the case of magnetically driven outflows as treated by 
Lovelace, Berk, and Contopoulos (1991, hereafter LBC). 
In the LBC model the outflows come predominantly from 
an annular inner region of the disk of radius ~ Ra where 
the disk rotation rate is lj$ = y 'GM / 'R A . Thus we as- 
sume that the outflows come from a region of the disk 
which is approximately in Keplcrian rotation. For the 
present situation, shown in Figure 1, it is clear that we 
must have n < Ra < fi ■ Further, we will assume Ra 
is close in value to r\ with (Ra — 1. For condi- 

tions where the outflow from the disk is relatively low tem- 
perature (sound speed much less than Keplerian speed), 
equations (16) and (18) of LBC imply the general relation 
E out = u)q ~ 3GM M out I (2Ra) ■ This equation can be 
used to eliminate L out in favor of E out in equation (6). 
Recalling that M out = M accr we have 



(l-S^)E out + 



3GMM a 
2^T 



-(1-Si 



(7) 



where S = Ra/^2 < 1- 

The energy dissipation in the region of the disk r = r\ 
to r 2 heats the disk and this heat energy is transformed 
into outgoing radiation E rad . Thus we have 



E. 



rad 



r2 , GMM dlsk {r) 
dr - 



(8) 



(Shakura 1973; Shakura & Sunyacv 1973), where 
M disk (r 2 ) = Maccr and M disk {ri) w 0. The essential 
change in M d i sk (r) occurs in the vicinity of Ra so that 
E rad w {3GMM accr /2)(R A 1 - r^ 1 ). Thus equation (7) 
becomes 



(l-6*)E out = (w*-W2)T 1 - 



3GMM a 
2R7 



1(1 -Js). (9) 



Thus the power from the spin-down of the star w*Ti must 
be larger than a certain value in order to drive the outflow. 

For magnetically driven outflows, the value of E ou t can 
be written as 



En 



^ ui R\ B<f MSc 



(10) 



(equation 34 of LBC), where T§ is a dimcnsionlcss numeri- 
cal constant ~ 0.234, and Bo is the poloidal magnetic field 
at the base of the outflow at r — Ra- We take the simple 
estimate Bq = n/R\ , which omits corrections for example 
for compression of the star's field by the inflowing plasma. 

Next we consider the torque on the star T\. Because 
most of the matter inflowing in the accretion disk at 
r = r 2 is driven off in outflows at distances r > n, 
the stress is necessarily due to the magnetic field. The 
magnetic field in the vicinity of r\ has an essential time- 
dependence owing to the continual processes of stellar flux 
leaking outward into the disk, the resulting field loops be- 
ing inflated by the differential rotation (LRBK), and the 



reconnection between the open disk field and the closed 
stellar field loops. The time scale of these processes is 

t\ ~ 2irri/vK(ri). We make the estimate of the torque 
Ti = -2wrl (2Az)(B r B (j> )i/(ATr), where Az is the vertical 
half-thickness of the region where the magnetic stress is 
significant, and where the angular brackets denote a time 
average of the field quantites at r ~ r\ . The magnetic field 
components, B ri _B^, with B^ oc —B r necessarily, must be 
of magnitude less than or of the order of the dipolc field 
B\ = [ijr\ at r = r\. The fact that B$ has the opposite 
sign to that of B r is due to the fact that the B^ field arises 
differential rotation between the region r < n, which ro- 
tates at rate w*, and the region r > ri, which rotates 
at rate ujk{ti) < u>#- Also, it is reasonable to assume 

Az ~ n. Therefore, 



2 _ - A 1 



Ti w a. x r\B{ = a 



R? A 



(11) 



where a < 1 (the time average of a(t)) is a dimension- 
less constant analogous to the a— parameter of Shakura 
(1973) and Shakura & Sunyaev (1973). (Note that be- 
cause n ~ Ra, ot and a\ = (ri/i?^) 3 ^ are of the same 
order of magnitude.) 

Substituting equations (10) and (11) into equation (9) 
gives 



a ri 



R 2 a 



A 'c r 



Ra 



5^ 



;(1 



Here, we have introduced two characteristic radii 
corotation radius, 



(12) 



the 



(GM\ 



1/3 



1/3 p 2/3 



1.5 x 10 s cm M^'°P : 



with Pi = (27r/o)»)/ls the pulsar period and Mi = 
M/MQ. [For a young stellar object, r cr « 1.36 x 

10 12 My 3 P 1 2 ( ^, where P\c, d is the period in units of 10 
days.] The second is the nominal Alfven radius 



r A = 



M accr VGM 



2/7 



3.6 x 10 8 cm 



4/7 
^30 



where the accretion rate Mn = M accr /(10 17 g/s) with 
10 17 g/s w 1.6 x lQ- 9 M Q /yr and ^ 30 = ^/(10 30 Gcm 2 ) 
with the magnetic field at the star's equatorial surface 
n/r 3 = 10 12 G(10 6 cm/r) 3 . [For a young stellar object 

r A « 1.81 x 10 12 cm^3/ 7 5 /(M 2 / 7 M 1 1/7 ), where the normal- 
ization corresponds to a stellar radius of 10 11 cm, a sur- 
face magnetic field of 3 x 10 3 G, and an accretion rate of 
1.6 x 10 _8 M Q /yr.] The corotation radius is the distance 
from the star where the centrifugual force on a particle 
corotating with the star (w 2 r) balances the gravitational 
attraction (GM/r 2 ). The Alfven radius is the distance 
from a non-rotating star where the free-fall of a quasi- 
spherical accretion flow is stopped, which occurs (approxi- 
mately) where the kinetic energy-density of the flow equals 
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Fig. 2. — The plot shows the effective Alfven radius Ra (normalized by the nominal Alfven radius ta) as a function of the 
normalized corotation radius r cr /rA for a — 0.2 obtained from equation (12). We have neglected the 5 3 / 2 = (Ra/^) 3 ^ 2 
terms compared with unity and taken the value !Fq = 0.234 from LBC. Only the part of the solid curve above the dashed 
line Ra = r cr is consistent with our assumptions. The dashed horizontal line r to indicates the "turnover radius" of the 
disk rotation curve or effective Alfven radius in the regime discussed by LRBK where M accr falls onto the star and the 
star spins-up. For this line the turbulent magnetic diffusivity of the disk is taken to be aD = 0.1. The points abed 
and the associated vertical lines represent possible transitions between spin-down and spin-up of the pulsar as discussed 
in the text. 



the energy-density of the star's dipole field. Note that the 
assumptions leading to equation (12) require Ra > r cr . 

Notice that Ra (or n ~ Ra) is the effective Alfven ra- 
dius for a rotating star. It depends on both ta and r cr in 
contrast with the common notion that the Alfven radius is 
given by r a even for a rotating star. ^From equation (11), 
the spin-down rate of the star is I(duj*/dt) = — a fx 2 /R 3 A , 
where I is the moment of inertia of the star (assumed con- 
stant). Thus the spin-down rate depends on both ta and 

Figure 2 shows the dependence of Ra on r a and r cr 
for a sample case. For conditions of a newly formed disk 
around a young pulsar, the initial system point would be 
on the upper left-hand part of the curve. Due to the 
pulsar slowing down (assuming fi and M constant), the 
system point would move downward and to the right as 
indicated by the arrow. In this region of the diagram, 

Ra ~ \fc T t r 7 J 4 /rcl 4 <x y/uj (for S -C 1), so that the torque 

on the star is T x w (GMM) 3 / 2 /(/i\/5 ^ 3/2 )- Thus the 
braking index is n = —3/2, where n is defined by the re- 
lation = —const w™. Numerically, 

1 dP _ 1.55 x IP' 5 P^ /2 (M 1 Mnf 12 ^ 
P dt yr y/E fx 30 £45 

where J45 = I/(10 45 g cm 2 ). For Ra r cr , the mass ac- 
cretion rate to the star is small compared with M accr , but 
some accretion may occur due to 'leakage' of relatively 
low angular momentum plasma across field lines near r\ 
(Arons & Lea 1976). 

As Ra decreases, the spin-down torque on the star 
increases. Over a long interval, Ra will decrease to a 
value larger than r cr but not much larger. In this limit, 
mass accretion to the star M\ may become significant. 
Our treatment can be extended to this limit by not- 
ing that Mi = M 2 - M out , L x = MiUJir 2 - T u and 
Ei = Mi(—GM/2ri)—u*Ti, where T\ is given by equation 
(11). In this limit the accretion luminosity is GMMi/r*, 
where r* is the star's radius. Figure 2 is not changed ap- 
preciably for Mi < M accr . 



Further spin-down of the star will cause the system 
point in Figure 2 to approach the right-most part of the 
curve. Further spin-down of the star is impossible. At 
this point of the evolution, the only possibility is a tran- 
sition to the spin-up regime. In this regime the effective 
Alfven radius is the 'turnover radius' r to of the disk rota- 
tion curve calculated by LRBK, the star spins-up at the 
rate Idui^/dt = M accr [GMr to ) x ^' 1 ', and most of the disk 
accretion M accr falls onto the star. The dashed horizontal 
line in Figure 2 indicates rto which is necessarily less than 

r cr • 

The location of the turnover line r to in Figure 2 suggests 
the possible evolution shown by the sequence of points 
a-^fe— >c-^d— > a. The system can jump down from 
point a where the star spins-down to point b where it spins- 
up. The spin-down torque at a is —otfi 2 /R Aa (where Racl 
is the effective Alfven radius at point a) , whereas the spin- 
up torque at b is M accr [GMr to Y^ 2 ■ The magnitude of the 
ratio of these torques is 

spin-down = JlA^ l% fRAa\ 1/2 (u) 
spin -up \RAaJ \ Tu, J 

With the system on the r to line, it evolves to the left. Be- 
cause r to < r cr , there must be an upward jump from point 
c to point d. For this case the torque ratio is given by 
equation (14) with i?^ a — > RacL- From point d the system 
evolves to the right. For the example shown in Figure 2, 
the torque ratio is w 3.45 for a — > b whereas it is w 0.87 for 
c — > d. The vertical line c — > d is at the left-most position 
allowed for the considered conditions, but this transistion 
could also occur if the line is shift to the right. The line 
a — > b can be displaced slightly to the right or it can be dis- 
placed to the left to be coincident with the c — * d line. In 
the latter case the torque ratio for the spin-down to spin- 
up jump is approximately equal to the torque ratio for the 
spin-up to spin-down jump and is w 0.87. The smaller the 
horizontal separation of the a — > b and the c — > d lines, 
the shorter is the time interval between jumps. 

Summarizing, we can say that the horizontal locations 
of the transistions, a — > 6, and c — > d, are indeterminate 



MAGNETIC PROPELLER OUTFLOWS 



5 



within a definite range. The locations of the jumps in the 
(Rai r cr) plane may in fact be a stochastic or chaotic in 
nature and give rise to chaotic hysteresis in the of the spin- 
down/spin-up behavior of the pulsar. The jumps could 
be triggered by small variations in the accretion flow (M 
for example) and magnetic field configuration (the time- 
dependence of a in the torque T\). Analysis of the accret- 
ing neutron star system Her X-l (Voges, Atmanspacher, & 
Scheingraber 1987; Morrill et al. 1989) suggests that the 
intensity variations are described by a low dimensional de- 
terministic chaotic model. The transistions between spin- 
down and spin-up and the reverse transistions may be de- 
scribed by an analogous model. 

The allowed values in Figure 2 have r cr /rA < const = k, 
where k ~ a 2 / 3 . This corresponds to pulsar periods 

For some long period pulsars such as GX 1+4 this in- 
equality points to magnetic moment values ^30 apprecia- 
bly larger than unity. Periods much longer than allowed 
by (15) can result for pulsars which accrete from a stellar 
wind (Bisnovatyi-Kogan 1991). [For a young stellar object, 

equation (15) gives P < 8d(a/M 1 5/7 )(^ 6 5 /Mi 8 ) 3 / 7 .] 

3. DISCUSSION 

This work presents a new investigation of the 'propeller' 
regime of disk accretion to a rapidly rotating magnetized 
star. The work considers the field configuration proposed 
by LRBK, the theory of LBC on magnetically driven jets, 
and the conservation of mass, angular momentum, and 



energy to derive an expression for the effective Alfven ra- 
dius Ra (equation 12) and the spin-down torque on the 
star Ti (equation 11). Our work is in qualitative accord 
with that of Li and Wickramasinghe (1997) who also con- 
sider the propeller effect of Illarionov and Sunyaev (1975). 
Our Figure 1 is similar to Figure 4 of Li and Wickramas- 
inghe for the spin-down regime, and our earlier work on 
the spin-up regime (LRBK) agrees with their Figure 3. We 
find that Ra depends not only on /i, M, and M, but also 
on the star's rotation rate W*. Because Ra decreases as 
decreases, there is a minimum value of uj* or a maximum 
value of the pulsar period P = 2tt/oj*. The model points to 
a mechanism for 'jumps' between spin-down and spin-up 
evolution (and the reverse transition). In our picture, in 
a spin-down to spin-up transition, for example, the effec- 
tive Alfven radius decreases by an appreciable factor going 
from > r cr to < r cr . The propeller goes from being "on" to 
being "off" in this transistion, which is a change between 
two possible equilibrium configurations. The transistions 
may be stochastic or chaotic in nature with triggering due 
to small variations in the accretion flow or in the magnetic 
field configuration. The ratio of the spin-down to spin-up 
torques (or the ratio for the reverse transition) is found to 
be of order unity (equation 14). This agrees with obser- 
vations of for example Cen X-3 (Chakrabarty et al. 1993) 
and GX 1+4 (Chakrabarty et al. 1997; Cui 1997). 
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